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We are going to analyze the nucleation of vortices in a thin mesoscopic superconducting disk. 
The Gibbs free energy from London model is formulated. This energy function, with an arbitrary 
configuration of vortices, is assoicated with the disk's size, applied magnetic field and finite tem- 
perature. Then, the optimal solution is obtained by differentiating with respect position r, for 
fixed applied field and temperature. We also investigate the stability of the different vortex states 
inside the disk and compare our results with those of other theoretical studies and with available 
experimental observations. Our results agree with experiments. Besides, we formulate the Gibbs 
free energy by Ginzburg-Landau (GL) model. Theoretically, the free energy by GL model takes the 
superconducting density into account but not in London model. Our simulations from both theories 
show the same quantum states of vortex. We find that the Gibbs free energy by GL model is smaller 
than by London model. 
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I. INTRODUCTION 

It is an interesting topic that how vortex nucleates in- 
side mesocopic supcrfiuids and superconductors for few 
decades. In our studies, we are going to study the vor- 
tices on the small disk. We formulate the Gibbs free 
energy in terms of the applied magnetic field H, the 
size of sample and temperature. The London model 
and Ginzburg-Landau model are applied in our stud- 
ies. The formulation of Gibbs free energies are combined 
with the kinetic energy, magnetic energy and supercon- 
ducting density energy. We also modify the Gibbs free 
energies by adding the term of entropy. It means that 
the order-disorder of vortices inside the sample has been 
taken into account. We found that this Gibbs free energy 
by Ginzburg-Landau model can be aprroximated into the 
one by London model. Of course, the London theory only 
describes the macroscopic phenomena, while the GL con- 
siders the microscopic influences at the same time. After 
that, we will show that the GL model and London model 
agree with the experiments. 

It is better to have an overview of some backgrounds 
of superconductors in mesoscopic scale. The nucleation 
of vortex in superfluid films and small superconductors 
as a function of flow velocity (superfluids) and applied 
magnetic field (superconductors) have been the subject 
in recent years. It is well known that the characteristics 
of small systems [l[ are affected by the sizes and the ge- 
ometries. As the size of the system becomes small, the 
boundary conditions become more important. In addi- 
tion, the vortex interactions in mesoscopic scale samples 
are different from those in large samples. Geim and co- 
workers found that there could be paramagnetic 
Meissner effect in small superconductors. Chibotaru et 
aZ.[5( and Mel'nikov [6| studied the possibility of anti- 
vortices and multi-quanta vortices penetrating thin meso- 
scopic square samples. Schweigert et a/.@, Q concluded 



theoretically that the multi- vortex state would transform 
into a single giant vortex state as the magnetic field and 
the disk thickness increase. On the other hand, Okayasu 
and co-workers |9| reported that the giant vortex could 
not be explained according to their experiments. 

Most recent theoretical studies can explain 

the creation of different vortex states and their stabil- 
ity well. In our studies, we modify the current results 
of small thin disk and compare with experimental obser- 
vations [161 ] . Experimentally, Grigorieva and co-workers 
[l6j studied vortex configurations in mesoscopic super- 
conducting disks, and found that the vortices form con- 
centric rings or shells, rather like shell filling in atoms and 
nuclei. They found that the configrations of vortices (the 
so-called magic numbers), corresponding to the appear- 
ance of new shells, are reproducible in many experiments 
for different applied fields H and diameter D of their 
disks. There are some theories which explains the mech- 
anism for vortex shell filling. However, we would like to 
modify the formulation of free energy by introducing the 
idea of entropy. Theories [HJ also predict stable configu- 
rations well generally. Nonetheless, in some cases, they 
are different from those observed experiments (l6| . It is 
appropriate to remark here that while experiments are 
performed at finite T ^ OK temperature, most theoreti- 
cal studies are only valid at T = OK. 

Finally, we'll extend the preliminary work of Sobnack 
et al. [l0| to include temperature by taking into account 
of the entropy energy associated with the vortex state. 
Using GL model and London model [l7l - l24l |. we write 
down the free energies of the disk with an arbitrary con- 
figuration of vortices arranged in shells. The free energy 
is minimized and the optimal vortex configurations are 
then obtained. We also found that Gibbs free energy in 
GL model is less than from London model. 
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II. FORMULATION AND METHODOLOGY 

A. Configuration of Mesoscopic Disk 

Assume that there are vortices inside a thin disk under 
th e appli ed magnetic field H with radius R and thickness 
d (l9l424l | (volume is V — ttR 2 cI). The small disk in an 
applied magnetic field is under H = H\s. = V x A app , 
perpendicular to the plane of the disk. We restrict the 
scale to the case i?<A = A 2 /d(Ais usual London pene- 
tration depth), and d <C A, with H near the lower critical 
field H c i. We further assume that d is in the magnitude 
of r c (where r c is the radius of the vortex core), so that 
the sample can be approximated to 2D plane. We then 
choose cylindrical coordinates (r, 8, z) to simplify our cal- 
culations and we apply the method from Buzdin [ToL [TTj 
and Sobnack [9j. The restriction on R implies that the 
screening effects of the superconducting currents j s are 
suppressed: the whole disk can be thought of as soaking 
in the applied field H, so that the local magnetic field 
B = V x A is approximately the same inside and out- 
side the sample, B « H. We formulate the problem at 
a temperature T < T c which is in superconducting state, 
where T c is the critical temperature. 



B. Ginzburg-Landau Model 

In this part, we will study the Gibbs Free energy [lH, 
|2~H in our small disk. We assume that the electrons in 
the superconducting state are in couple (cooper) pairs. 
From the I s * quantization, the order-disorder of these 
electron pairs can be descibed as wave function [l9l - l2"lj ]. 
^(x) = I -0 1 e* 6 * , where refers to the norm and 9 is 
the phase of the wave function. We can think this as a 
vector with planar orientation. We also assume that the 
spatial variation of the \ip\ to be very small. Under zero 
applied magnetic field, the free energy density [H, [H| of 
the system becomes 



9 = go + 



higher order terms, (1) 



where a and (3 are parameters associated with supercon- 
ducting material, and go corresponds to the reference free 
energy density. This is the form of Talyor's approxima- 
tion. If we further consider that the applied magnetic 
field H which is added on the system, two more signif- 
icant terms need to concern. The first one is similar to 
the Hamiltonian of cooper pair under the electromag- 
netic field where the classical dynamic equation of 
the electrons pair is 



d 2 r 

m ^2 =? E +; vxB - 



whereas E is the electric field which is zero in our case; 
v is the velocity of superconducting current and q = 2e. 



Indeed, this is the kinetic enegry of electrons under the 
magnetic field. In the view of quantum mechanics, 
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A represents the magnetic vector potential and we can 
choose an arbitary gauage that is invariance. The sec- 
ond term is about the magnetic energy stored inside the 
system, (B — H) 2 , where we define B as a local field 
generated by the superconducting system itslef due to 
the Meissner effect. This magnetic term refers to the 
normal energy in the sample. The total Gibbs free en- 
ergy density under applied magnetic field H, according 
to the Ginzburg-Landau theory, is hence, 
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(2) 



For a small disk, the screening is weak. The term of 
magnetic energy is small comparable to the kinetic en- 
ergy and H « B so that this magnetic term can be ne- 
glectable later. 

Next, we need to minimize the given Gibbs free energy 
in order to find out the optimal solution. By variation 
principle, derivative with respect to the order parame- 
ter we obtain the first equation of Ginzburg-Landau 
theory: 



a* + ^*|*| 2 + -!- ( -V - -A 
2m \ 1 c 



* = 0, (3) 



with boundary condition (f V — |A) - n = that normal 
component equals to zero, and only tangential compo- 
nent remains jl5l |26|. The 1 st equation is analogy to the 
non-linear schrodinger equation which is time indepen- 
dent. While we take the differentiation with respect to 
A to the free energy, the 2 nd GL equation comes out. It 
corresponds to supperconducting density current j s : 

j s = — (**W - *V**) - — ***A, (4) 
2mi mc 

which is similar to the quantum density current flow. 



C. Vortices inside the Small Thin Disk 

Up to this point, we are going to solve out the solu- 
tion of vortices inside the thin disk. Suppose the order 
parameter is ty(r) = I'S/le 1 ' 5 , as mentioned above. Then, 
we consider the gradient of order parameter iKS/^ which 
is associated with the quantum momentum, and related 
to the kinetic energy component. In the cylindrical co- 
ordinate, 
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By the way, the disk can be treated as very thin so that 
we can approximate and then map it into 2-D cylindrical 
coordinate. We assume that the amplitude of \^>\ is am- 
lost constant except at some singularities and the edge 
of the disk. It is on the other side to assume the spatial 
variation of |W| is very small and hence, 



1 r)v[> „ j 

V* w -^770 = -W*0. 

r oO r 



(5) 



Putting back into the Gibbs free energy equation, we 
obtain: 



0, 
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(6) 



What we need to do next is to handle the kinetic energy 
term. The gradient of the phase k = V$ represents the 
number of circulations of the vortex, and thus k is the 
winding interger number. 

Up to this point, we need to further discuss the physi- 
cal mechanism of vortices in small type-II superconduct- 
ing sample. If the applied field H is smaller than the 
first critical field H c i, the sample will tend to oppose H 
(the Meissner effect). In the range of H c \ < H < H C 2, 
applied flux lines will penetrate into the disk and then 
vortices are being nucleated, with each vortex carrying 
a flux <p v = l<j> , where / G Z (I > for vortices and 
I < for anti-vortices) and 4>o — he/ q is the quantized 
fluxion. Typically, a core of radius r c in a vortex is in 
nano- magnitude of order 10 _9 m. On the other hand, 
outside the region of core radius, it more belongs super- 
conducting state. Within the core region, it is a normal 
state rather than a superconducting state. Two more 
pararmeters are needed to clarify the range of normal 
and superconducting states: London penetration depth 
A and coherent length £. A measures how the distance 
relates to the applied magnetic field down to neglectable 
small inside the superconductor. And £ measures how the 
range of sample from normal state, affected by applied 
field, back to superconducting state. We can imagine 
a vortex in the water pool. There is the circulation flow 
around the center which is usually a small hole. The hole 
can be regarded as the normal state while the flow refers 
to the super-current. In the thin mesoscopic disk, the 
circulating superconducting current around each vortex 
is approximated as } v oc -6. 

Let's think about the kinetic energy which is con- 
tributed by the superconducting current, 
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It is obvious that A v = (f>ok/2Trr is the magnetic po- 
tential generating by vortex. Now the Gibbs free energy 
expression is hence, 

9-90 = ,>|*| 2 + ||*|* 

After this, we have to handle the first two Talyor's terms 
of the free energy: 



11 2 1 1 2m ^ e 2^\a 



where we define £ 2 = h 2 /2m\a\ [H, For the sake of 
minimization, I'J'ol 2 = M//3, then the Taylor's terms will 
become: 
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where the dimensionless parameter is = 1 i'/| , i , o|- By 
intcrgating both parts J d 3 r and we define free eenrgy 
from GL mdoel, 
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(9) 



where G = J gd 3 r. Besides, we can also approximate the 
applied magnetic potential A app to be very close to local 
magnetic potential A [IFJ, [ll| . 

Similar the treatment to previous result [T3, [ll| , un- 
der certain temperature and magnetic field, the Gibbs 
free energy can be re-written in terms of number of vor- 
tices inside the small thin disk. We define L = — 1, 0, 1 as 
a vortex at the center and N as a number of vortices on a 
shell of the disk. If there is no vortices, say L — 0, TV = 0, 
the normalized order parameter is |#| oc c, where c is a 
constant. Suppose the effect of the edge is neglectable. 
This parameter describes how much the superconduct- 
ing and normal regions inside the system. While the 
square of order parameter represents the superconduct- 
ing density of the system, n s oc \^\ 2 cx d (c' is constant). 
We only need to take into account where there are some 
vortices which will cause the superconducting region be- 
comes normal state. Hence, we could assume that region 
(volume) of normal state created by a vortex can be ap- 
proximated as 7r£ 2 d, where £ refers to the coherent length 
of a vortex. The order parameter can be thus assumed 



|*| 2 cx C '(l-(£/i?) 2 (A + L)), (10) 
where N -\- L is the total number of vorticities. 
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D. Vortex Shell Configurations on the Disk 

The first two terms of the free energy are related to 
the density of superconducting sample. So what we have 
to concern now is the kinetic energy term, 



h 2 \* \ 

2m 

ft 2 |*o| 2 
2m 



— If A, 



(A v - A 



appy 



(11) 



Here we assume the super-density n s oc |^| 2 varies 
slowly in spatial and A ~ A app . From this point, 
we can assume the local field be B ~ H = (0,0, H), 
whereas the magnetic potential vector can be chosen as 
A = (—yH/2,xH/2,0) - Landau gauage. More specifi- 
cally, 



A = -B x r = - 

2 2 



i j k 
H 

x y z 



Now we need to choose the optimal postion of vortices. 
It is obtained by minimization of total free energy with 
respect to r. If there is only one vortex, the optimal 
position will be situated at the center of the disk. If 
there is more than one vortices, the vortices will at the 
position circulating around the center with radius r = fj. 
The main concern is the interaction between vortices and 
local magnetic field. We can approximate this as analogy 
with electrostatics, by adding each vortex of flux fa at 
ri in the disk, then there is an image anti-vortex of flux 
— fa at r' ; = (R/ri) 2 ri beyond the disk and leads to 



A v = ^[$ i (r-r i )-$,(r-r' i )]0, 



(12) 



with $i(r) = (f)i/(2nr). In our formulation, the vortices 
can be filled up to two shells. We assume that N\ and 
N 2 be the number of vortices occupy on the two shells; 
where L and M are the number of quanta of vortices at 
the center and in the shell respectively. 

After some calculations, we can show that, at T = K, 
the Gibbs free energy of a configuration of N\ vortices at 
r, = ri and N 2 vortices at rj = r 2 , each of flux <j) = M(j>o 
in the shell, and a concentric vortex (anti- vortex) of flux 
<j> = Lfa (L > for a vortex and L < for an anti- vortex) 
as 



2m 



S 2 |* | 2 27rd 



g GL {L,N u N 2 ) 



"|pH*| 2 + i|l*! 4 + 9L ond (L,N 1 ,N 2 )\*\ 2 , (13) 
whereas 

gLand^m,^) = -h 2 + L 2 \n--2LN l M\nz 1 

4 r c 

- 2LN 2 Mlnz 2 -Lh + g'(N 1 ,0) 
+ g'(N 2 ,0)+g 12 (N 1 ,N 2 ), (14) 



where g = mG/(h 2 \^>Q\ 2 ird) is the dimensionless Gibbs 
free energy with reference, h = HirR 2 j fa is the dimen- 
sionless applied field and z% = Ti/R (i — 1,2) shows 
the optimal position. Furthermore, we discover that the 
dLond is in fact the Gibbs free energy from London model. 
g' Lon d{Ni, 0) and g' Lond (N 2 , 0) are the dimensionless free 
energies of Ni and N 2 off-centre vortices respectively, 
with 

g'(Ni, 0) = NiM 2 In — - N^N, - l)M 2 In z t 
+N,M 2 ln(l - z 2 ) - N % Mh(l - z 2 ) 

^r h liM^W (15) 

2 ^ 4sin 2 (™/iV l ) 

for (i = 1, 2) and, as is usual, we have introduced the core 
radius r c as a cutoff whenever ri = rj. The interaction 
energy gi 2 (Ni, N 2 ) of the N\ vortices in the first shell 
(radius r\ ) and the N 2 vortices in the second shell (radius 
r 2 ) is 

g 12 {N u N 2 ) = 
2 ^ 1 + z\z\ - 2z\z 2 cos[a + 2n(n/Ni - m/N 2 )\ 
q n z \ + z\ - 2z 1 z 2 cos[a + 2ir(n/N 1 - m/N 2 )} ' 



(16) 



where n £ [1, Ni—1], m S [1, ^2 — 1] and a is the misalign- 
ment angle between vortices in the two shells. Typically, 
we will take L = M = 1 for the reason of mininization of 
energy. 



E. Entropy of Ordered Vortices associated with 
Temperature 

Next, we consider the arrangement of vortice inside the 
mesoscopic thin disk. The entropy [lpj is introduced in 
our formulation. It measures how the order and disorder 
of the vortices placing on the disk. As the disk is in 
mesoscopic scale, the system is not allowed to put too 
many vortices on it. In our case, we assume there are two 
shells (rings) on the disk. The physical interpretation is 
that there are a number of arrangements of vortices on 
the shells, so that entropy exists when temperature is 
greater than zero but small than the critical limit. The 
general idea is thus 



G(L, N!,N 2 ,T) = G(L, N U N 2 )- TS, 



(17) 



where S = ks ln(VF) is total entropy of vortices and ac- 
cording to the Boltzmann statistics, W is the number 
of configurations. At finite temperatures T ^ OK, one 
has to take into account the entropy S associated with 
N\ vortices at = t\ and N 2 vortices at = r 2 . This 
gives an additional term — ksT{h\Wi + lnW^), where 
Wi — 2irr;j2Nir c (i = 1,2), giving the dimensionless free 
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energy as 
9gl(L, Ni,N 2 , t) 



t = 



= g GL (L,Ni,N 2 )-t(2hnr + 2\n R/r c 
- lniVi -\nN 2 +lnzi + lnz 2 ), (18) 



where zi = ri/i?, z 2 = r 2 /R and i = inksT / {h 2 \$> \ 2 itd) 
is the dimensionless temperature. The magnetization M 
of the disk follows from d{G GL +M -H)/dH = and this 
gives the reduced magnetization m as m(L, Ni, N 2 ) = 
-dg GL (L,N 1: N 2 ,t)/dh. 



III. SIMULATION RESULTS AND ANALYSIS 
A. Discusion of cut-off Dimensionless Parameters 

We will go through the cut-off value of dimensionless 
parameters briefly. The normalized parameters t, £ and 
A, and their effective ranges will be focused on. For 
the London equation, the penetration depth can be ex- 
pressed as A = y/mc 2 j A-n:n s q 2 , [24[ ■ The normalized 
temperature is t = mfcsT / (h 2 n s ird), where n s — \^o\ 2 - 
Let's consider the parament from GL formaulation, while 
l/irn s = 4\ 2 q 2 /mc 2 , t = A 2 fc B f67r 2 T/0Qfi, whereas (f) = 
hc/q. This exactly matches with what we find from the 
London theory. The effective A e // oc X 2 /d w 0.6 ~ Q.9R 
[l9| is chosen as the weak screening effect by the small 
disk. It can be understood that magnetic field lines pass 
through and occupy in most regions of mesoscopic scale. 
The dimensionless ratio n = A/£ is an indicator to de- 
fine type I and II superconductors. Here as we make use 
of effective penetration depth, we will assume the effec- 
tive coherent length in parallel, such that £ e ff oc ^\/d. 
The physical interptation of £ e / / should be linking to the 
normal area state, Tr^ 2 ff , created by a single vortex. 



B. 



Vortices inside Disk by Ginzburg-Landau 
Theory with Entropy 



For a given temperature i, we minimize 
gch{L, Ni, N 2l t) with respect to z\ and z 2 for a 
range of applied magnetic fields h, with different L, Ni 
and N 2 . Then the optimal Gibbs free energy is obtained 
and the transition for each state is considered. 

The critical temperature of Niobium is around T c = 
9.1 K, so that in our dimensionless units, the critical 
temperature is t c — 0.7; while t = 0.14 corresponds to 
the operative temperature T — 1.8 K at which the ex- 
periments were performed. Figures 1 and 2 show Gibbs 
free energy for t = and t = 0.14 respectively. The 
parameters d, R, A are chosen so that our disk corre- 
sponds to the sample studied experimentally by Grig- 
orieva et al. fl6l ] and some approximations. [Nb, radius 
R w 1.5 ~ 2.5 urn, A w 90 nm, £ w 15 nm, d w 0.1 ~ 0.3£] 

As h increases from zero (Figure 1), the free energy of 
the screening currents increases quadratically with h as 
expected until the first critical field hi ~ 5.3 is reached 
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FIG. 1: The free energy g with respect to h at t = is shown. 
The stable configurations are successsively (0,0) — > (1,0) — ► 
(0,2) (0,3) -> (0,4) -»■ (0,5) -> (1,5) (1,6) (1,7) -> 
(1,8) ->• (0,2,8). The states (L,JV) = (1,5), (1,6), (1,7), 
(1,8) and (0,2,8) (corresponding to total flux 6<jf>o, 7cio, 8<fto, 
9<f>o and 10<f)o respectively) are more stable than (N, L) = 
(0,6), (0,7), (0,8), (0,9) and (0,10) or (1,9). 



when a single vortex [(L,N) = (1,0)] penetrates the 
disk at the centre (see Fig. 1). This state persists un- 
til the second critical field h 2 ~ 7.4 is reached at which 
the single center-vortex splits into two off-center vortices 
[(0,2)]. As h is further increased, more off-center vor- 
tices nucleate on the first ring, forming successively a 
triangle, a square and a pentagon, until the sixth crit- 
ical field he ~ 12.9 is reached. Then it is energetically 
more favorable for the next vortex to nucleate at the 
center of the disk [(1,5)] than to form a hexagon of six 
off-center vortices with a vortex at each vertex. This re- 
sult agrees with those of other studies 0, 0, H HH, G3 
and is also analogous to the result of the study by Yarm- 
chuck el al. [14| on the nucleation of vortices in super- 
fluid He, which showed that a central vortex would ap- 
pear in the system. As h increases further, further off- 
center vortices enter the disk and nucleate on the first 
ring, the stable vortex states going through transitions 
(1,5) (1,6) (1,7) -> (1,8), until the tenth critical 
field ~ 17.0 is reached when, instead of the tenth vor- 
tex nucleating to form state (1,9), the vortices rearrange 
themselves to form the state (L, Ni,N 2 ) — (0, 2, 8), with 
no vortex at the center of the disk, two vortices on the 
first ring and eight on the second. The entry of each ad- 
ditional vortex at successive critical fields is accompanied 
by a jump in the magnetization m of the disc (shown in 
Figure 2). 

At t = 0.14, the Meissner state persists until the ap- 
plied magnetic field reaches the first critical field hi ~ 5.3 
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FIG. 2: The magnetization m with respect to h at t = 0. The 
stable configurations are (0,0) — > (1,0) —5- (0,2) — > (0,3) — > 
(0, 4) -> (0, 5) -»• (1, 5) -► (1, 6) -»• (1, 7) (1, 8) -► (0, 2, 8). 



t = 0.14 
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FIG. 3: The free energy g as a function of h at t = 0.14 
(T = 1.8 K). The stable vortex states and transitions are 
(0,0) -> (1,0) -+ (0,2) (0,3) -> (0,4) -> (0,5) -> (1,5) -+ 
(1,6) -¥ (1,7) -)■ (0,2,7) -> (0,2,8) as h increases (for the 
range of h shown). 



when a single vortex [(L,N) = (1,0)] nucleates at the 
center of the disk (see Figure 3). At the second critical 
field hi ~ 6.7, the energetically favorable configuration 
(with total flux 2</>o) is the state with two off-center sin- 
gle vortices (L, N) = (0, 2). As ft. increases further, more 
vortices penetrate the disk, with the fluxiod state going 



FIG. 4: The magnetization m as a function of h at t = 0.14 
(T = 1.8 K). The stable configurations are (0,0) -> (1,0) -s> 
(0,2) -> (0,3) -> (0,4) -» (0,5) -»• (1,5) -> (1,6) -> (1,7) -» 
(0,2,7) -»• (0,2,8). 



successively through the transitions (0, 2) — > (0, 3) — > 
(0,4) -> (0,5) (1,5) (1,6) -)■ (1,7) until the next 
critical field h$ w 14.9 is reached. Then an extra vor- 
tex enters the disk, but the stable vortex state with total 
flux 90o is the state (0,2,7): no vortices at the centre 
of the disk, and the nine vortices form two rings, with 
two vortices on the inner ring and seven on the outer. 
As h is further increased, a further vortex penetrates the 
disk and nucleates on the outer ring, forming the state 
(0,2,8). These results are in direct agreement with the 
experimental observations of Grigorieva et al. [16( . Fig- 
ure 4 shows the magnetization of the disk in which dis- 
continuous jump is found. 

Figure 5 extracts parts of figure 1 respectively, giving 
that Gibbs free energy with total flux 9cf>o. The solid 
line and triangle in Fig. 5 give the free energy of the 
vortex state (0,2,7) and (0,3,6) correspondingly (two 
shells configuration); while the star represents the state 
(1,8), totally with flux 9^o- The most stable state at 
t = (T = OK) (Fig. 5 left) is the state (1,8), with 
a central vortex and eight vortices on a shell. The next 
stable comfiguration is (0,2,7) and the unstable one is 
(0,3,6). 

Whereas at t = 0.14 (T — 1.8 K) (Fig. 5 right), the 
most stable state is the state (0, 2, 7), with no vortex at 
the center of the disk, two vortices on the inner ring and 
seven on the outer ring. This result is in very good agree- 
ment with the experiments of Grigorieva et al. [l6| (who 
found that at T — 1.8 K, the state (1,8) was observed in 
only just a few cases, while the state (0, 2, 7) was, by far, 
the most frequently observed state) . From the studies of 
Baleus et al. [l5[, they theoretically predicted only the 
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FIG. 5: Free energy with total magnetic flux = 94>o- Blue, 
red and green lines represent the states (1,8), (0,2,7) and 
(0, 3, 6) correspondingly. Left diagram reveals curves at t = 
while right diagram is at t = 0.14. 
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FIG. 6: Total flux = 90o: State (1,8) with one concentric 
vortex and eight off-center vortices. 
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FIG. 7: Total flux = 9^ : State (0,3,6) with three off center 
vortices in the 1 st ring and six off-center vortices in 2 nd ring. 
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FIG. 8: Free energy with total magnetic flux = 1000 • Blue, 
red and green lines represent the states (1,9), (0,2,8) and 
(0, 3, 7) correspondingly. Left diagram reveals curves at t = 
while right diagram is at t — 0.14. 



state (1,8) will be present. Figs 6-7 show the supercon- 
ducting density with states (1,8), (0,2,7) and (0,3,6) 
sequently. The color bar, representing the the norm of 
super- density, is shown from interval [0,1]. We can see 
that there are some holes inside the disk. Physically, 
those holes mean that the superconducting states are be- 
ing destroyed by the field lines, which turns to normal 
state with zero magnitude. Magnitude with one regards 



as the superconducting region without destroy. 

Figure 8 shows the corresponding result for vortex 
states with total flux lO^o- The triangle shows the free 
energy that depends on h with state (0, 3, 7), the star that 
represents the state of (1,9) and the line represented by 
red colour is (0, 2, 8). At t = (Fig. 8 left), the most sta- 
ble vortex state is the state (0, 2, 8), with states (0, 3, 7) 
and (1,9) having almost the same (slightly higher) en- 
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FIG. 9: Total flux = lO^o: State (0,2,8) with two off center 
vortices in the 1 st ring and seven off-center vortices in 2 n 
ring. 




FIG. 10: Total flux = 10^ : State (0, 3, 7) with two off center 
vortices in the 1 st ring and seven off-center vortices in 2 nd 
ring. 



ergy. 

At t = 0.14 (T = 1.8 K) (Fig. 8 right), the state 
(0, 2, 8) is again the most stable state, with the other two 
states having higher energies. This result is in very good 
agreement with the experimental studies of Grigorieva 
et al. [IH wn0 reported that the state (0,2,8) was the 
most frequently observed state and that the state (1,9) 
was never observed in their experiments. Baleus et al. 



15], on the other hand, predicted only the state (1,9) 
for the state with total flux 10(fo- Figs 9 and 10 are the 
density in states (0, 2, 8) and (0, 3, 7) respectively. 

C. Thermal Fluctuation of Quantized Vortices 

Experiment ally (lq . it is found that the state (0,2,7) 
with total flux 10</>o occurs more probable than the state 
(1,8). The similar case appears in the subsequent state 
(0,2,8) and (0,3,7). Statistically, the occurance of vor- 
tex state (0, 2, 7) is around 7 times more than (1, 8); while 
the appearance of (0, 2, 8) is about 2-3 more than (0, 3, 7). 
In order to explain the situation, we would like to re- 
ansonably approximate the probability of vortex state 
by the Maxwell-Boltzmann distribution. It is to say, 
p oc exp^—Gch/kBT), as the probability of certain con- 
figuration, where Gi is the Gibbs free energy at the state 
i and ks is the Boltzmann constant. For dimensionless, 
the probability expression can be written as 

p oc exp(-g GL /t), (19) 

where g G L is the normalized free energy. Our approxima- 
tion shows that the occurance of state (0, 2, 7) is around 
10 times more than (1,8). We also predict that state 
(0, 3, 6) will appear 2 ~ 3 times more than (1, 8), in which 
it does not show in the experiment. Besides, the appear- 
ance of state (0, 2, 8) is around 2 ~ 3 times more than 
(0,3,7). Both cases are consistent in the experiement. 

D. Gibbs Free Energy Compared with London 

Approximation 

Finally, we compare the results between the Ginzburg- 
Landau and London theories. The Gibbs free energy de- 
rived from GL model can be approximated to the one 
by London model. From the general results, we found 
that both Gibbs free energies produce the same vortex 
state. One can find that the free energy from London 
model (figllip has a higher level of magnitude than the 
one by GL model. It is due to the fact that the GL model 
also takes the supperconducting density |^| 2 = n s into 
account. As the temperature or magnetic field increases 
further, it will undoubtedly destory part of the super- 
density which turns to normal state. Th free energy has 
therefore got a fall. 



IV. CONCLUDING REMARKS 

We have formulated the Gibbs free energy by 
Ginzburg-Landau equations in a mesoscopic disk. A crit- 
ical review of the early study of Sobnack and Kusmartsev 
[l(J, which the London equations are applied, is made. 
We also introduce the idea of entropy into the Gibbs 
free energy. Inclusion of the temperature term —TS (by 
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FIG. 11: The free energy gLond as a function of h at t = 0.14 
(T = 1.8 K) .The stable configurations are (0,0) ->■ (1,0) -> 
(0,2) -»• (0,3) (0,4) -> (0,5) -► (1,5) -► (1,6) -» (1,7) 
(0,2,7) -> (0,2,8). 



taking into account the entropy 5 associated with the 
non-center vortices) lowers the free energy of some of 
the vortex states and stabilizes them. Our results are in 
agreement with those of the recent experiments of Grig- 
orieva and co-workers (lfjj . It is found that our results 
matches with those of Baleus et al. 15] in many cases, 
and there are only some states, (1,8) and (1,9) which 
show some differences. Those different states are also 
predicted by our modified formulation which coincide in 
the experiments [l6| . A possible reason for the disagree- 
ment between the theory of Baleus et al. [To] ] and ours 
is that the experiments are performed at finite temper- 
atures r ^ 0K, whereas the study of Baleus et al. [H[ 
concern the circumstance at T = K. The free energies 
by London and Ginzburg-Landau equations have been 
investigated. We find that free energy formulated by 
Ginzburg-Landau theory provides a lower level of energy. 
Physically, it means that the ingredient of superconduct- 
ing density is taken into account. More flux lines will 
break the order symmetry and some regions will transit 
to normal. Broadly speaking, both free energies can be 
obtained the same configuration of vortices in small disck. 
It is shown that the free energy by Ginzburg-Landau the- 
ory can be reasonably approximated to London theory. 
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